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Abstract 

The concept of time-coarsened density matrix for open systems has 
frequently featured in equilibrium and non-equilibrium statistical me¬ 
chanics, without being probed as to the detailed consequences of the 
time averaging procedure. In this work we introduce and prove the 
need for a selective and non-uniform time-sampling, whose form de¬ 
pends on the properties (whether thermalized or not) of the bath. It is 
also applicable when an open microscopic sub-system is coupled to an¬ 
other finite system. By use of a time-periodic minimal coupling model 
between these two systems, we present detailed quantitative conse¬ 
quences of time coarsening, which include initial state independence of 
equilibration, deviations from long term averages, their environment 
size dependence and the approach to classicality, as measured by a 
Leggett-Garg type inequality. An interacting multiple qubit model af¬ 
fords comparison between the time integrating procedure and the more 
conventional environment tracing method. 

PACS number(s): 03.65.Yz,75.10.Jm, 05.30.-d 


1 Introduction 

Though not usually presented as such, the density operator formalism pro¬ 
vides at least a partial resolution of the quantum mechanical time-arrow 
conundrum (time inversion invariant equations versus the expected and ob¬ 
served uni-directional trend towards equilibration), and does so by invoking 
some averaging procedure that spoils the unitary development of the pure 
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state density matrix (DM). The equilibration issue has been particularly 
acute for closed systems (the Universe or part of it hermetically segregated 
from the rest, such as trapped ultracold atoms) and has featured in many 
articles both much cited ones m-m and those quite recent m-m- The com¬ 
monality among these appears to be the inclusion of a random process that 
is not implicit in the defining equation. There is no consensus on the depen¬ 
dence (or on its absence) on initial conditions of the equilibrating system 
and in what situation this (or that) occurs le.g..|12jl. 

The time-arrow conundrum lends impetus to the present work which de¬ 
rives density matrices in arbitrary representation for open systems, such that 
a microscopic system is coupled to an ancilla (which can be another micro¬ 
scopic system or a finite or infinite ’’bath”), and does so by a method of time 
averaging of the microscopic system’s state function, in contrast to other, 
more frequently employed procedures, in which a stochastic process or trac¬ 
ing over the ancilla’s states is invoked m- It is a development arising from a 
previous work in which justifications of the time averaging method were 
given based on the ergodic hypothesis (equating time- and state sampling- 
averages) |15( fTB] , and the ” minimal coupling model” (shown in equation ([T]) 
below) was heuristically derived. [We note that this model, also known as 
the ” monochromatic coupling model”, has quite a history starting with laser 
coupling m and other topics, supported by detailed theoretical analysis of 
the model in the context of Floquet’s theorem by m and m, and very 
recently in relation to Landau-Majorana-Stiickelberg-Zener interferometry 
[20], also [HI) 

After describing (for completeness and in a nutshell) our simplified em¬ 
ployment of the time averaged DM method (in section 2), the present work 
makes the following advances requisite for the establishment and practice of 
the method: It shows first (in section 3) that the method satishes the basic 
theoretical conditions on DM. It then identifies quantitatively the forgetting 
of initial conditions (section 4). In the short, but central section 5 (propped 
up by the formal proof in Appendix D) we delineate, apparently for the first 
time, how the ’’time coarsening” method for the DM is to be used for open 
systems coupled to finite-sized environments and, by a natnral extension 
to infinite-sized environments, giving rise in general to non-uniform time 
weighting. Continuing, the paper develops an interacting N-qubit model (in 
section 6.1) in order to show how the standard deviation (’’distance”) of the 
open system DM vanishes with the size of the ancilla. This leads us (in 
section 6.2) to a quantitative comparison between the ancilla tracing and 
time averaging procednres regarding the ancilla’s size dependence of the DM 
standard deviation in each procedure. The formalism is then applied in sec- 
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tion 6 to a Leggett-Garg type inequality to investigate quantitatively how 
the inequality, violated by quantum systems, is reinstated as the ancilla size 
approaches macroscopics. To make the paper self contained we have briefly 
reiterated in Appendixes A and B the bare essentials of [H]. As already 
noted, the way time averaging is to be done for an open (microscopic) sys¬ 
tem coupled to a finite size environment is the subject of the (mathematical) 
Appendix D, with implications also for general macroscopic environments. 

2 Time Averaging in the Minimal Coupling Model 
(MCM) 

This model was introduced in with motivations and heuristic justifica¬ 
tion given there, in the form of the Hamiltonian 

H{t) = Eaz + k(TxS\u{ijjt) {oj —>■ 1) (1) 

and the (transposed) solution of the relevant time-dependent Schrodinger 
equation with pre-fixed initial conditions 

Wf' = ( 2 ) 

are Pauli matrices, 2E is the Zeeman splitting, k the coupling strength 
of the spin to the reservoir and the time [t) dependent circular function rep¬ 
resents the minimal effect of the latter on the former; u, I label the upper and 
lower component in the spinor solution. The DM pij in an ij representation 
is derived from the time averages over a time window 2At of 

Mij{t)=<i\fi{t) ><'tjj{t)\j > (3) 

in the form 

-[ rt+At 

^ ^ Jt-At ^ 

3 Necessary Conditions for Density Matrices pij 

These were formulated by Fano, as follows [22| : 

1. Matrix is Hermitian. This is clearly the case in the above MCM. 

2. Matrix trace is one. For the time averages this follows from the nor¬ 
malization to unity of the wave function at any time. 
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3. In any representation the diagonal elements pu > 0: this follows from 
the definition. 

4. When diagonalized by a unitary transformation Trp^ < 1. This re¬ 
quires that for each diagonalized eigenvalue {DM)n{p) of the DM, 
0 < {DM)n{p) < 1, in which the equalities hold for pure states. In the 
present spin-half formalism with two DM eigenvalues this translates 
to 

0<i(l±^^)<l (5) 

where the unity in the parentheses comes from unit trace. Here, with 
a little manipulation, 

1 /*£-t-At 

G{t) ^ [— 

2At Jt-M 

1 ^t+At 

+ ( 6 ) 

At Ji-At 

(involving the upper (u) and lower {1) components of the wave func¬ 
tion) and this has to be numerically less than 1 for the last Fano re¬ 
quirement to hold. Now, G{t) is clearly non-negative and is maximal 
when the two components are throughout real and positive, and limits 
for this quantity need to be sought. The derivation of a maximum of 
unity is given in Appendix C. 

This completes the proof for the satisfaction of the Fano requirements 
for time averaged DM. 

4 Equilibration 

The first question that needs to be asked about the model is how it de¬ 
scribes equilibration characteristics of the spin system. (The mode and 
time-duration of the evolution to the equilibrated state, discussed recently 
in e.g. [23], is outside the scope of the present time-integrated approach. His¬ 
torically, the dependence of the rate of equilibration on the power-spectrum 
of the coupling seems to have been first given in [24|.) Defining the equili¬ 
brated value of Pij{t) as that value which is only minimally dependent on 
time (an issue discussed further in the ’’Distance” section of this work), we 
ask how does this depend on the initial conditions and on the strength of the 
coupling kl Also, at what values of this coupling and of other parameters is 
’’full equilibration” (= the maximal entropy stage) attained? Figures 2 and 
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Figure 1: Upper diagonal density matrix (DM) as function of the coupling 
strength k for the half-Zeeman splitting E = 1, and the following starting 
values: 1, 0.75, 0.5. 



Figure 2: Upper diagonal DM as function of the coupling strength k 
for a different half-Zeeman splitting E = \/2, and the following starting 
values: 1, 0.75, 0.5. 
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3 show the results. The curves show that the starting value is forgotten 
beyond A: = 2 in figure 1 with E = \ and beyond about A: = 3 in figure 2 
with E = \/2- 

Another important lesson from the figure is the strong fluctuations in 
the results, so that for physical interpretation only the rough averages and 
tendencies are meaningful. This strongly fluctuating property of the MCM 
and its analogues are well known from previous works m-mm- 

5 Non-Uniform Time-Weighting 

5.1 Finite Bath: Discrete Energy Levels 

For a finite, though large, environment (bath) its level scheme is strictly 
discrete and the density of states is singular. The coupled subsystem’s in¬ 
stantaneous energy E(t) is continuous. The two systems form together a 
pure state with eigenenergy E^otai ([S]). Then the continuous (Lebesgue) 
time integration in equation ([5|) with uniform time weighting has to be 
replaced by a discrete sum. Alternatively, the integral has to weighted by 
the singular factor 

Y^S{Etotal - E{t) - Ei)/ < 'ljj{t)\6{Etotal - Hsystem “ Ei)\'lp{t) > (7) 

i 

where the sum runs over all energy values of the environment and 5{x) is 
the Dirac delta function. A proof of this formula is given in Appendix D 
and an application is provided in section 6. 

5.2 Macroscopic Bath: (Quasi-) Continuous Energy Density 

For this case, in a straightforward extension of the previous expression, the 
weighting factor for time coarsening is 

D{Etotal - E{t))/ < 'll:{t)\D{Etotal - Hsystem)\'4^{t) > ( 8 ) 

in which D{Etath) is the distribution of bath energies. This form was used 
in there based on heuristic reasoning, for a bath distributed as in a 
canonical ensemble. 

6 ’’Distance” from equilibrium state 

Pursuing the environment-tracing formalism to derive the reduced density 
matrix of a small subsystem, and starting with an arbitrary but representa¬ 
tive pure state of the subsystem plus the environment, in their above quoted 
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paper the Bristol group have shown that the long term time averaged de¬ 
viation (’’Distance”) of the subsystem state (its truncated density matrix) 
from the equilibrated state is (or is upper-limited) proportional to the square 
root of the ratio between the subsystem’s dimension in its Hilbert space and 
that of environment (Equation 8 in [U], also 0 ). Before establishing an 
analogous result for the time integration formalism we show that an inter¬ 
acting N half-spin model verifies numerically the predicted inverse square 
root relationship for the ’’distance” [26]. 

6.1 iV-spin model 

Such systems with finite Hilbert state dimensions have been extensively 
studied for entanglement and other properties. A recently built processor 
consisting of eight interacting qubits and subject to time dependent inter¬ 
action was investigated as a means to achieve feasible quantum computing 

m- 

In our model the total system of N 1/2-spins (qubits), consisting of one 
subsystem (labelled 0) and (N — 1) bath 1/2-spins {i = 1...N — 1), interact 
according to a Hamiltonian 

H(A^) — ^ zi 1 “ ^ ) Qij^xi^xj (9) 

with energy parameters, so chosen as to avoid any degeneracies in the 
eigenenergies, and two sets of coupling parameters for gij = gji, given by 

Eo = l;Ej = 1.5sin2^; 5oi=2.5,3.5; =-75,1.5; l<i,j<N-l 

( 10 ) 

The matrix of this Hamiltonian, that has to be diagonalized to obtain its 
eigenenergies and pure-state eigenfunctions, is of a square dimension of 2^ 
X 2^. To formally write out in the spin’s representation the matrix for large 
N and then let it be solved by a routine (e.g., Fortran or Mathematica), we 
have used the following trick, not used (to our knowledge) heretofore: We 
have replaced the ordinal, decimal number of the matrix rows (or columns) 
in the routine by its binary representation and have formalized the matrix 
elements in this representation. Thus, there is an entry at the decimal 
(13,25), or (25,13) matrix position with the value of 53^5 = 55 ^ 3 , because 
in the binary representation we have for N = 8 (with the spins labelled as 
0,1,..., 7, in which the last 7 digits label the bath spins) 

13 = [OOOOlIOl], 25 = [OOOllOOl] (11) 
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Figure 3: Eigenenergies (in arbitrary units) against the serial number r of 
increasing energies in an interacting A^-spin system for N = 5,6,7 


with the 3 and 5 spins that flip and counter-flip identified by a superbar. 

The 2^ ordered energy values (labelled with the ordinal number index r) 
obtained by numerical diagonalization of the Hamiltonian equation ([9]) are 
shown in Figure [3l The figure shows sparseness in the lower range, uniform 
linear increase in the middle and some super linear increase towards the 
upper end. 

The reduced density matrix of the small subsystem (of size 2 x 2) is 
obtained in the bath-tracing procedure by tracing over 2^“^ bath states the 
total system’s density matrix for any of its pure eigenstates (whose number 
is 2'^). The ’’distance” is taken as the root mean square deviation of p^ over 
the different choices from the 2^“^ bath eigenstates. 

6.1.1 N-spin distances 

In Figure m we show our computed ’’Distances”, plotted logarithmically 
(Logio) against the dimension of the Hilbert space of the ( — 1) half¬ 

spins with which the system (also: a half-spin) interacts. The full curve, 
the long broken line curve, the medium broken line curve and the short 
broken line curves are for different parameter sets in the Hamiltonian; the 
one (reaching up to A^ = 9) is for a parameter set in which all the axC^x 
coupling strengths have the same value 2.5 for every (bath-spin) - (system- 
spin) coupling goi and the value of 0.7 for all inter- bath couplings pij. The 
other curves are for varied coupling strength of the same order. The energy 
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LogDistance 



Figure 4: ’’Distances” versus N, explained in the text. 


scale is set by the Zeeman-splitting (= 2) for the subsystem; the bath spins 
have varying splittings of similar magnitudes. 

The two straight, dotted lines bordering the computed curves show pu¬ 
tative dependencies of the ’’Distance” on the Hilbert space dimension of 
2 ^'^“^), the lower dotted line following (decimal logarithmically) the inverse 
of this dimension and the upper dotted line the inverse square root. Asymp¬ 
totically, the computed values appear to follow the inverse square root law, 
in line with the predicted upper limit dependence in Eq. ( 8 ) of [ 6 j. 

6.2 DM through time summation with random partitioning 

As noted in section 5 and Appendix D, for a finite-sized ancilla the time 
integration over a window of 2At needs to be replaced by a discrete time- 
summation. We have subdivided the window into P segments (partitions) 
having randomly chosen relative lengths and summed the DM eigenvalues at 
P points corresponding to the mid points of the segments. Figure [5] shows a 
characteristic set of results [28]. For finite P the discrete summation corre¬ 
sponds to a finite bath while for P —)• oo the Standard Deviation (SD) should 
approach the time integration value, corresponding to an infinite bath envi¬ 
ronment. The computed curves appear to show initially a P “2 relationship, 
similar to the inverse square root dependence for the environment tracing 
result in figure 01 followed for large P by one like P~^. 

The two figures, 0] and O show the congruity of the environment trac¬ 
ing and the time integration methods for a general result in |6| which is 
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AvDistance 



Figure 5: Distance (or Standard Deviation SD of the density matrix upper 
eigenvalue over a large time interval) against the number P of randomly 
partitioning the time integration interval 2At = dvr and summing the value 
of the upper eigenvalue of the density matrix at the midpoint of each par¬ 
tition. The curves are for parameter values of the Hamiltonian shown in 
equation ([1]) oj = l,e = 0.5. The solid lines are for /c ~ 1 and the broken 
lines k = 1.5, with two different random partitioning in each case. The two 
dotted straight lines are, in the logarithmic plotting of the figure, slopes of 
and P~^, respectively. 
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model independent. 


7 A Leggett-Garg-Wrachtrup (LGW) Inequality 

The Bell inequalities for correlations between measurements at different lo¬ 
cations, famously potentially violated by quantum systems while observed 
by classical systems, were formulated for temporally subsequent measure¬ 
ments by Leggett and Garg [29]. A review of experimental and theoretical 
developments is found in m- A modification which tests the system for 
’’reality” ( meaning the characterization of the physical system, irrespective 
of whether it is measured or not) and ’’stability” (’’the conditional probabil¬ 
ity Pr{i, t|f, to) to find the system at times t, to in the same state i, provided 
that this probability depends only on the time difference”) was proposed in 
m by an inequality written, with to = 0 being the starting time and t a 
subsequent instance, as 

Pr(i, 2t|t, 0) — Pr^(t, t|i, 0) > 0 (12) 

We use this inequality to answer two questions: 

1. Does the (reduced) density matrix arising from an open microscopic 
quantum system interacting with an ancilla have the property of a 
quantal or a classical system, in the sense of violating or observing 
the above inequality? (We recall that in a measurement of a quantal 
observable, after interaction with the measuring apparatus, the trun¬ 
cated DM has the status of classical probabilities.) 

2. Assuming a classical status for the DM of an open system linked to 
an ’’infinite” environment and a quantum status for the same system 
when linked (entangled) to another microscopic system, is there a way 
to decide between the status of the combined systems in terms of the 
requirements of the above LGW inequality? 

It is proposed that the time integrating formalism leads to answers in the 
following sense: The execution of the time integration is dependent on the 
coupled ancilla. For a macroscopic ancilla this has to be carried out over 
an effectively full temperature range (or over a full period, if there is a 
periodicity in the system), whereas for a finite ancilla the integration is to 
be restricted, either to a summation or to a limited range of the integration. 
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We then test the violation of the LGW inequality by increasing the time 
averaging range 2At from zero (the pure state case) to its effectively full 
averaging value 2At = Att, appropriate to a macroscopic bath. 

We solve equation ([2|) with the initial condition = 0) = 1, V'z(O) = 0. 
Then the inequality in equation (|1‘2I) , takes the simple form 

|V’n(2t)|2 _ > 0 (^3) 

The results for some choices of the parameters in equation ([1]) are shown 
in Figure [ 6 ] for an infinitesimal averaging window, 2At —)• 0, corresponding 
to an instantaneous or pure state. It is seen that for a small range of 
time values [when the left hand side (LHS) of equation (fT3]l takes negative 
values] the inequality is violated. Such violation may either indicate the 
quantal nature of the system or to the breakdown of the pre-conditions in 
m- We next select some time values where each curve is most negative 
and average the probabilities around these values, continuously increasing 
the extent of the integration window. As seen in Figure [71 above some 
finite size time windows the LGW inequality becomes satisfied, though it 
was violated for an infinitesimal time window (the pure state case). The 
transition from negative to positive values may indicate the transition from 
quantal to classical nature of the integrated DM in a continuous manner, 
though this is also contingent to the satisfaction of the ” stability” hypothesis 
in m- For the present purposes the results demonstrate the facility with 
which application of the time averaging procedure straddles the quantal- 
classical gap. 

8 Conclusion 

While the main result of this paper is the derivation of non-uniform weight¬ 
ing for the time coarsening process in section 6 , a discussion of its application 
is appropriate. It may be argued that the Minimal Coupling Model MCM 
used in this work is too specific to represent the wide variety of real life 
interactions between a micro-system and its environment. (Extensions of 
the model, possible tasks for future work, are described in Appendix B.) 
However, it is almost axiomatic that, excluding the small proportion of in- 
tegrable systems, in real life (exemplihed in [ 6 ] by the Cooling Coffee Cup), 
thermal or equilibration processes depend only in a minor way on the de¬ 
tails of the micro-system-environment coupling, and indeed most modelings 
of this coupling have hinged on their amenability to solution, rather than 
on their being a true description of real life situations. The strength of the 
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InstantProb 



Figure 6: Pure sate results for the LHS of equation (fT^ obtained from 
solutions of the time dependent Schrodinger equation with a choice of pa¬ 
rameters, each choice drawn with different lines. The curves are mainly 
positive, but there are some times at which they are negative, thereby vio¬ 
lating the inequality in equation (fT^ 
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Figure 7: Time averages of probabilities as functions of the time averaging 
window 2At. The time choices are those, at which the curves in the previous 
figures are most negative, with the drawn curves corresponding to those in 
the preceding figure. 


time-averaging approach in this paper is the ease with which it yields results 
(resolution of the time arrow problem being perhaps the most obvious one); 
it also opens the way to explore concrete physical situations by the selective, 
non-uniform time averaging, as set out in section 5. 

9 Acknowledgements 

Thanks are due to Rafael Ruppin for substantial help in the calculations 
and to Ronnie Kosloff for a discussion. 


A Original outline of the three steps to construct 
the density matrix 

In the quoted previous work [TTj it was shown how the decohered-truncated 
phase of the full subsystem-environment state (their DM), leading into the 
reduced subsystem state, can be obtained by a time averaging approach. 
The abstract steps needed to achieve this and the main assumptions behind 
it were given by [14] in an ’’Outline”, which is briefly reproduced in this 
Appendix. The Hamiltonians for which the system’s states were determined 
were (i) a Rabi-model (a single-spin in interaction with a single vibrator). 
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(ii) two mutually coupled spins, each coupled to a single vibrator. For com¬ 
pleteness these Hamiltonians, formulated in a semi-classical, time-dependent 
language, are also shown in this article in Appendix B. While it is clear that 
a Hamiltonian in which a time dependent term appears represents an open 
system, what our results have shown is that even with the arguably simplest 
form of the time dependent term (namely, a single subsystem-environment 
harmonic interaction term), one obtains results equivalent to those in the 
standard environment-tracing formalism. The simplifying (and approxima¬ 
tive) steps necessary to pass from the Leggett et a/’s (spin-oscillatory ensem¬ 
ble) model [32] to ours were described in section 2.3 of the earlier article. 
The present text has extended the time-averaging formalism to issues not 
considered before, including second moments (fluctuations) of the DM (in 
section 6) and its approach through time averaging to classicality (section 

7). 

1. We adopt the von Neumann dehnition 

>< 'ipait)\j > (14) 

In this definition the summation index a represents the values of all 
coordinates, variables etc. external to the system (e.g., those of the 
environment affecting the system) and appearing also in the Hamilto¬ 
nian. Thus the set '4’a{t) for all a’s forms a time dependent ensemble 
of states. The variables of the system themselves are implicit (not 
written out) in V'a(i)- 

2. We solve only for a single external condition thus dispensing with the 
a index in the wave function, but obtain p{t) as the average over an 
adequate set of adjacent times: 

2 rt+At 

^ dr < i|V’(T) >< V’(T)|j > (15) 

This should be equivalent to equation (flTl) if the ergodic hypothesis 
holds for the duration 2At. While for several cases the off-diagonal ma¬ 
trix elements are small or vanishing, Pij{t), as defined above, represents 
in general a mixed state whose diagonalized form Pij{t) —>■ pf^{t)dij sat- 
ishes < 1 (section 5 in this paper). 
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B Minimal Coupling Hamiltonians 

(a) Single spin: 


H{t) = ecTz + kax sin(wt) 


(16) 



h{t) = '^[Eiazi + huzi cos{uit + a') + k[axi sm{ujit + a')] (18) 



having written in the first line the total Hamiltonian, comprising the parts 
in the next two lines. First, the spin energy terms in which Ei are energies of 
the spin systems, the cr’s are Pauli matrices operating in the respective spin 
spaces; k, a, together with their tagged partners, are parameters of the spin- 
boson couplings and ui is the frequency of the external source. This external, 
boson source is classical and for it the Hamiltonian need not be written out. 
The Hamiltonian does not include back-reaction on the source, which can at 
least approximately be justified for periodic coupling and energies. The last 
line is the spin-spin interaction term. A two-spin version of this interacting- 
spin Hamiltonian was treated in M- 

C Extremization of G{t) in equation ([6]) 

We now prove that the maximum of G{t) is unity when the wave function 
components ipu{t) and '^^(f) are constant during the 2At integration range 
m- By implication, G{t) is less than unity (mixed state case) when the 
components are genuinely time dependent. The method of proof is a calcu¬ 
lation of variations. 

It has been stated in text that a maximum of G{t) occurs when the wave 
function components are throughout real and positive. We further simplify 
by writing 



( 20 ) 


Written as 


1 /“tH-Af 1 /— 

G[f{t)] = {2 fit') - l)dt'f + [— fit')/il - fit'))dt'f 


( 21 ) 
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a variation gives 


5G[/(t)] = 


+ 



rt-\-At 1 rt-\-At 

/ / 4f{t')6f{t')dt'] 

Jt-At Jt-At 

rt-\-At j - 1 rt-\-At j - 

Jt-At ^ Ar Jt-At 


Pit') 

Vi-PiP 


)5fit')dt'] 


(22) 


Remarkably, the cofactor of Sf{t') vanishes when f{t) = /, a constant 
throughout the range of integration, so that G[f{t) = /] is then an ex¬ 
tremum. Its value is unity, independent of /. That it is also a maximum 
can be shown by evaluating G{t) for the case that the function /(t) takes 
two values /i ,/2 / fi in fractions F, (1 — F) (respectively) of the integra¬ 
tion interval and showing that the leading term in G[fi, f 2 ] is negative for 
I/ll, I/ 2 I < 1, namely 


^\ /i-/2)^[4 + (/i + /2)^][1 + 


4 + (/i + /2)2 

(4-(/i + /2)2)2 


(23) 


D Proof of Sum Formula in equation ([7j) for a 
Finite Bath. 


It is assumed that the subsystem and the large, but finite sized bath form 
together a microcanonical ensemble with mean energy of Etotai and small en¬ 
ergy uncertainty v. The range of the subsystem’s energies E{t), that vary in 
time, is divided up into N segments, numbered n(= 1,..., N) each of spread 
v. The segments contain Vn+i — Vn discrete (supposed non-degenerate) en¬ 
ergy levels of the bath, such that these bath state energy levels satisfy 

Etotai = E{t) + E^ (24) 

for some value of E{t)p E{tn)] situated within the segment spread zz. The 
corresponding time spread over the segment is tn+i—tn, having assumed that 
for the short passage time over the narrow segment spread the subsystem 
energy is a monotonic, single-valued function of time. This time spread is 
of the order of the bath’s relaxation time. 

The subsystem’s density operator (’’the reduced” density operator) p^it) 
is obtained in the bath tracing formalism by weighting the system’s proper 
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density operator pn{t) in each segment by the above number of bath states, 
giving 


Pn{t){rn +1 - Vn) 


/-E(Wi) 

Pn{t)Y. dE'5{Et,tal-E' -El) 

Tn 

rE+i dE(f'') 

E Pn{t) / dt'^:^5{Etotal - E{t') - El) 

rn ® 

r-n+l 

E] Pn{tr) (25) 

r=r„ 


having in the last line interchanged the order of the summation and the time 
integration and taken heed of the property of the Dirac delta function. We 
thus obtain for the weighted density operator in the time segment a discrete 
sum with the time sum going over all such times that the system energy 
is complemented to Etotai by a bath state’s energy. In the time-integrated 
formalism the system’s density matrix is obtained by integrating over time 
windows that are much wider than the time spread of a segment. Then 
the resultant matrix element is 


P^j ^ ^ ( Piji^r) (^ 6 ) 

where the time sum goes over all such times E within the time window that 
the instantaneous subsystem energy E{t) is complemented by a bath level to 
add up to Etotai- The proportionality constant is fixed, by the requirement 
that the trace of density matrix is unity. In the case of an inhnite bath with 
a continuous energy spectrum, the above sum is replaced by a time integral 
including the energy density of the bath levels D^[E^)\ 

pt+t-w 

Pij{t)(x dt'D\E^)p,j{t') (27) 

This formula is also applicable for the case that the bath is in thermal equi¬ 
librium at a temperature T, so that the bath energy density is proportional 
to 

e’^BT=e ^bt (28) 

by virtue of the expression in equation (|24l) for the bath energy. This result, 
originally due to [3l], has already been used in [14] , with the time indepen- 

dent factor e having been absorbed in the proportionality factor. 
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